Relativistic descriptions of final-state interactions in neutral-current neutrino-nucleus 

scattering at MiniBooNE kinematics. 
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The analysis of the recent neutral-current neutrino-nucleus scattering cross sections measured 
by the MiniBooNE Collaboration requires relativistic theoretical descriptions also accounting for 
the role of final state interactions. In this work we evaluate differential cross sections with the 
relativistic distorted-wave impulse-approximation and with the relativistic Green's function model 
to investigate the sensitivity to final state interactions. The role of the strange-quark content of the 
nucleon form factors is also discussed. 
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I. INTRODUCTION 

The MiniBooNE Collaboration has recently re- 
ported [l| a measurement of the flux-averaged differential 
cross section as a function of the four-momentum trans- 
ferred squared, = —q'^q^i, for neutral-current elas- 
tic (NCE) neutrino scattering on CH2 in a range up 
to « 1.65 (GeV/c)^. The neutrino-nucleus NCE reac- 
tion in MiniBooNE may be considered as scattering of 
an incident neutrino with a single nucleon bound in car- 
bon or free in hydrogen, but it can also be sensitive to 
contributions from collective nuclear effects, whose clear 
understanding is crucial for the analysis of ongoing and 
future neutrino oscillation measurements [H-Q . In addi- 
tion, NCE processes can be used to look for strange-quark 
contributions in the nucleon that may show up through 
the isoscalar weak current. 

Recent results on parity- violating electron scattering 
at Q2 = 0.1 (GeV/c)2 indicate that the strangeness 
contribution to the electric charge and magnetic moment 
of the nucleon is consistent with zero at 95% confidence 
level. In the axial form factor, under the usual assump- 
tion of a dipole Q^— dependence, the only free parameters 
within the relativistic Fermi gas (RFC) model 0, are 
the nucleon axial mass Ma and the strange-quark con- 
tribution As, determining the value of the axial form 
factor at = 0, that is related to the fraction of the nu- 
cleon spin carried by the strange quark. Recent chargcd- 
current quasielastic (CCQE) neutrino-nucleus measure- 
ments reported values of Ma « 1.2 - 1.3 GeV/c^, 
significantly larger than the world average value from the 
deuterium data of = 1.03 GeV/c^ [lllil- In agree- 
ment with these new results, the MiniBooNE NCE data 
provide a best fit for Ma = 1.39 ± 0.11 GeV/c^. 

A measurement of v{v)-'pvoion elastic scattering at the 
Brookhaven National Laboratory (BNL) [l3| at low 
suggested a non-zero value for the strange axial-vector 
form factor. However, the BNL data cannot provide us 
decisive conclusions when also strange vector form fac- 
tors are taken into account [31 ■ A determination of the 
strange form factors through a simultaneous analysis of 



vp, vp, and ep elastic scattering is performed in Ref. [15| . 

Since cross section measurements are a very hard ex- 
perimental task, ratios of cross sections have been pro- 
posed as alternative ways to search for strangeness 
Il7| . Moreover, they are expected to be less sensitive to 
distortion effects [3, [l^ ■ Taking advantage of the fact 
that at > 0.7 (GeV/c)^ single proton events can be 
satisfactorily separated from neutron and multiple nu- 
cleon events, the MiniBooNE Collaboration used the ra- 
tio of proton-to- nucleon (p/n) cross sections to extract 
the strangeness contribution to the axial form factor [l[ . 
The resulting value is As = 0.08 ± 0.26. Although af- 
fected by large systematic errors because of difficulties in 
the proton/neutron identification, this result is anyhow 
very interesting, since this is the first attempt to mea- 
sure As using the p/n ratio. In addition, exploiting its 
data sample of neutrino-nucleus events, the MiniBooNE 
Collaboration has also focused on the neutral-to-charged- 
current ratio [U, [3l , that can provide us useful and com- 
plementary information. 

The energy region considered in the MiniBooNE exper- 
iments, with neutrino energy up to « 3 GeV and average 
energy of « 0.8 GeV, requires the use of a relativistic 
model, where not only relativistic kinematics should be 
considered, but also nuclear dynamics and current opera- 
tors should be described within a relativistic framework. 
From the comparison with electron scattering data, the 
RFC turns out to be a too naive model to correctly ac- 
count for the nuclear dynamics, and thus the larger axial 
mass needed by the RFC could be considered as an ef- 
fective value to incorporate nuclear effects into the calcu- 
lation. Regardless from this question, a comparison be- 
tween the results of different models and the NCE Mini- 
BooNE data is important to clarify the role of the various 
ingredients entering the description of the reaction. 

At intermediate energy, quasielastic (QE) electron 
scattering calculations [2fl l21| , which were able to suc- 
cessfully describe a wide number of experimental data, 
can provide a useful tool to study neutrino-induced pro- 
cesses. However, a careful analysis of i^-nucleus NCE 
reactions introduces additional complications, as the fi- 
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nal neutrino cannot be measured in practice and a final 
hadron has to be detected: the corresponding cross sec- 
tions are therefore semi-inclusive in the hadronic sector 
and inclusive in the leptonic one. Several sophisticated 
models have been applied in recent years to i^-nucleus 
scattering reactions and some of them have been com- 
pared with the MiniBooNE data, both in the CCQE and 
in the NCE channels. At the level of the impulse ap- 
proximation (lA), models based on the use of a realistic 
spectral function [22l . [23| , which are built within a non- 
relativistic framework, underestimate the experimental 
CCQE and NCE cross sections unless is enlarged 
with respect to the world average value. The same re- 
sults are obtained by models based on the relativistic lA 
(RIA) [23 - [26| . However, the reaction may have signifi- 
cant contributions from effects beyond the lA in some 
kinematic regions where the neutrino flux for the ex- 
periment has si gnif icant strength. For instance, in the 
models of Refs. [27h31| the contribution of multinucleon 
excitations to CCQE scattering has been found sizable 
and able to bring the theory in agreement with the ex- 
perimental MiniBooNE cross sections without increasing 
the value of Ma- Fully relativistic microscopic calcula- 
tions of two-particle-two-hole (2p-2h) contributions are 
very involved and may be bound to model dependent as- 
sumptions. The part of the 2p-2h excitations which may 
be reached through two-body meson-exchange currents 
(MEC), in particular the contribution of the vector MEC 
in the 2p-2h sector, evaluated in the model of Ref. [3^ . 
has been incorporated in a phenomenological approach 
based on the superscaling behavior of electron scattering 
data [H,!!!. 

Within the QE kinematic domain, the treatment of 
the final-state interactions (FSI) between the ejected nu- 
cleon and the residual nucleus is an essential ingredient 
for the comparison with data. The relevance of FSI has 
been clearly stated in the case of exclusive (e, e'N) pro- 
cesses, where the use of complex optical potentials in 
the distorted wave impulse approximation (DWIA) is re- 
quired [13, [m, [35I - I41I . In the exclusive reaction, where 
the final state is completely determined, the absorptive 
imaginary part of the optical potential accounts for the 
flux lost to different final nuclear states. In contrast, 
the analysis of inclusive reactions needs all final-state 
channels to be retained, i.e., the flux must be conserved 
and the DWIA based on the use of an absorptive com- 
plex potential should be dismissed. Different approaches 
have been used to describe FSI in relativistic calcula- 
tions for the inclusive QE electron- and neutrino-nucleus 
scattering [i^ - lisj . Besides the relativistic plane- wave 
impulse approximation (RPWIA), where FSI are simply 
neglected, FSI have been included in DWIA calculations 
where the final nucleon state is evaluated with purely 
real potentials, either retaining only the real part of the 
relativistic energy- dependent optical potential (rROP), 
or using the same relativistic mean field potential con- 
sidered in describing the initial nucleon state (RMF). 
Although conserving the flux, the rROP is unsatisfac- 



tory from a theoretical point of view, since it relies on 
an energy-dependent potential, which reflects the differ- 
ent contribution of open inelastic channels for each en- 
ergy, and under such conditions dispersion relations dic- 
tate that the potential should have a nonzero imaginary 
term [i^. On the other hand, in the RMF model the 
same strong energy-independent real potential is used for 
both bound and scattering states. It fulfills the disper- 
sion relations [49| and also the continuity equation. 

In a different description of FSI relativistic Green's 
function (RGF) techniques [H-Iil HO, [HI are used. This 
formalism allows us to reconstruct the flux lost into 
nonelastic channels in the case of the inclusive response 
starting from the complex optical potential which de- 
scribes elastic nucleon-nucleus scattering data. Thus, it 
provides a consistent treatment of FSI in the exclusive 
and in the inclusive scattering and gives a good descrip- 
tion of (e, e') data [i^H^]. Moreover, due to the analitic- 
ity properties of the optical potential, the RGF model 
fulfills the Coulomb sum rule 0,111, [s^]. 

These different descriptions of FSI have been compared 
in [131 for the inclusive QE electron scattering, in [1^ 
for the CCQE neutrino scattering, and in [H^] with the 
CCQE MiniBooNE data, which are reasonably described 
by the RGF results without the need to increase the value 
of Ma and are generally underestimated by the other 
models. 

In this paper different relativistic descriptions of FSI 
for NC :^-nucleus reactions in the quasielastic region are 
presented and compared with the NCE MiniBooNE data. 
In these reactions a final nucleon is detected, like in the 
exclusive scattering, but since the final lepton cannot be 
detected, the final nuclear state is not completely de- 
termined and the cross section includes many possible 
final-state channels. In Refs. [l^ [H, [53l - [55| such a semi- 
inclusive scattering was treated with the same relativistic 
DWIA (RDWIA) model that was successfully applied to 
the exclusive (e, e'N) reaction, as a process where the 
cross section is obtained from the sum of all the inte- 
grated exclusive one-nucleon knockout channels. Results 
for both the semi-inclusive CC and NC cross sections 
were presented and the effects of possible strange-quark 
contributions on the cross sections and on other observ- 
ables were discussed. In RDWIA calculations the imagi- 
nary part of the optical potential produces a loss of fiux 
that accounts for the flux lost in each considered chan- 
nel towards other final channels. Some of these reaction 
channels are not included in the experimental cross sec- 
tion when an emitted nucleon is detected and for these 
channels this treatment of FSI is correct. There are, how- 
ever, other channels, which are excluded by the RDWIA 
approach but which can contribute to the semi-inclusive 
reaction. Some of these contributions may be small or 
negligible for the specific final state considered in the ex- 
clusive reaction, but may be much more important for 
all the final states of the semi-inclusive reaction. The 
fiux lost towards these channels can be recovered in the 
RGF, where the imaginary part of the optical potental 
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redistributes the strength in all the channels and the to- 
tal flux is conserved. The RGF, however, describes the 
inclusive process and, as such, can include also contri- 
butions of channels that are not included in the cross 
sections of the semi-inclusive reactions. Thus, according 
to the approach adopted to describe FSI, the RDWIA 
can produce cross sections smaller and the RGF larger 
than the experimental ones. The relevance of the contri- 
butions neglected in the RDWIA and added in the RGF 
depends on kinematics. It is not easy to disentangle the 
role of each specific contribution, in particular if we con- 
sider that both RDWIA and RGF calculations make use 
of phenomenological optical potentials, obtained through 
a fit of elastic proton-nucleus scattering data. 

In spite of these uncertainties, the comparison between 
the results of the RDWIA and RGF models with the 
MiniBooNE data can be helpful for a deeper understand- 
ing of the role of FSI in the analysis of NCE data. 



II. RESULTS AND DISCUSSION 

In all the calculations presented in this work the bound 
nuclcon states are taken as self-consistent Dirac-Hartrce 
solutions derived within a relativistic mean field approach 
using a Lagrangian containing a, uj, and p mesons [HgI - 
[ssj . The Energy-Dependent but A-Independent EDAI 
parameterization for the complex phenomenological po- 
tential of Refs. [G^I, which is fitted to elastic proton- 
^^C scattering data, has been used. The Energy- 
Dependent and A-Depcndent EDADl parameterization, 
which is fitted to clastic proton scattering data on sev- 
eral nuclei in an energy range up to 1040 MeV, has also 
been used for some calculations. We note that whereas 
EDADl is a global parameterization, EDAI is a single- 
nucleus parameterization, which is constructed to bet- 
ter reproduce the elastic proton-^^C phenomenology [60j . 
and also leads to CCQE results in better agreement with 
the MiniBooNE data 

In Fig. [T] the NCE (vN — > vN) cross sections aver- 
aged over the neutrino fiux arc shown as a function of 
and compared with the MiniBooNE data [l[. Cal- 
culations are performed in the RPWIA and RDWIA 
with Ma = 1.03 and 1.39 GeV/c^ and neglect possible 
strange-quark effects. A larger value of Ma gives a larger 
cross section, because of the dominant role played by the 
axial-vector current in NC scattering, but the enhance- 
ment is not linearly proportional to Ma and, therefore, 
also the shape of the cross section is slightly modified. 
In the comparison with data, the RPWIA results show 
a generally satisfactory, although not perfect, agreement 
with the magnitude, while some differences are obtained 
with respect the shape of the experimental cross section. 
In the RPWIA, however, FSI are completely neglected. 
The RDWIA results are smaller than the RPWIA ones 
and therefore also smaller than the experimental data. 
This is due to the imaginary part of the optical poten- 
tial, that in the RDWIA gives an absorption that reduces 
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FIG. 1. NCE flux-averaged {vN — Y uN) cross section as 
a function of calculated with two different values of the 
axial mass in the RDWIA (solid and dashed lines) and RP- 
WIA (dotted and dot-dashed lines). The data are from Mini- 
BooNE [H. 



the calculated cross section. In Fig. [T] the RDWIA cal- 
culations with the EDAI potential generally underesti- 
mate the NCE cross section, unless > 0.8 (GeV/c)^ 
for Ma = 1.03 GeV/c^ and > 0.6 (GeV/c)^ for 
Ma = 1.39 GeV/c^ We have checked that close RD- 
WIA results are obtained with the EDADl potential, but 
for small values of , where there are visible differences 
that can be attributed to the different imaginary parts 
of the two optical potentials in the low energy region. 

In Fig.[2]we show our RDWIA results with As = -0.18 
and -HO. 34 and Ma = 1.39 GcV/c^. These are the upper 
and lower limits for As found by MinibooNE [l[. The 
rcsuhs for As = with Ma = 1.03 and 1.39 GeV/c^ 
are also shown again for a comparison. The MiniBooNE 
NCE cross section is nearly independent of As, as the 
combined effects on proton and neutron events almost 
cancel. In the calculations a negative As produces an 
enhancement of the proton and a suppression of the neu- 
tron contribution, which arc approximately of the same 
size (see also [II])- In the case of positive As the effect is 
reversed. As a consequence, the effects of different values 
of As = arc minimal and smaller than the uncertain- 
ties due to Ma, whose value has a visible impact on the 
MiniBooNE NCE cross section. 

In Fig.|3]wc show our RGF results calculated with both 
EDAI and EDADl potentials and compare them with the 
RDWIA and the rROP results calculated with the EDAI 
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FIG. 2. NCE flux-averaged {uN — i- i/N) cross section as a 
function of Q^, calculated in the RDWIA with three different 
values of As and Ma = 1-39 GeV/c^. Results with As = 
and Ma = 1-03 GeV/c^ are also shown. The data are from 
MiniBooNE 



FIG. 3. NCE flux-averaged {i/N — > i>N) cross section as a 
function of calculated with the RGF-EDADl (solid hne) 
and RGF-EDAI (dashed line). The dotted and dot-dashed 
lines are rROP and RDWIA results calculated with the ED AT 
potential, respectively. The data are from MiniBooNE 



potential. All these calculations have been performed 
with Ma = 1.03 GeV/c^ and As = 0. The RGF cross 
sections with both optical potentials are larger than the 
RDWIA and the rROP ones. The rROP result, where 
a purely real optical potential is used, underestimates 
the experimental cross section for < 0.6 (GeV/c)^. 
We observe that a rROP calculation with a larger axial 
mass, e.g., AIa ~ 1.3 — 1.4 GeV/c^, is able to reproduce 
the data with good accuracy. However, we note that, 
independently of its result in comparison with data, the 
rROP model, which is based on an energy-dependent po- 
tential, has important physical drawbacks. The RDWIA 
cross section with the EDAI potential is the same already 
presented in Fig. [1] and gives in Fig. [3] the lowest cross 
section. The differences between the RGF results calcu- 
lated with the two optical potentials are clearly visible, 
although not too large, the RGF-EDAI cross section be- 
ing in good agreement with the shape and the magnitude 
of the experimental cross section, and the RGF-EDADl 
below the data only at the smallest values of con- 
sidered in the figure. The differences between the RGF 
results obtained with the two phenomenological optical 
potentials can give an idea of how uncertainties in the 
determination of this important ingredient can affect the 
predictions of the model. These differences are basically 
due to the different values of the imaginary parts of the 
two potentials, particularly for the energies considered 



in kinematics with the lowest values of Q^. The real 
terms of the optical potentials are very similar for differ- 
ent parameterizations and give very similar results. In 
the rROP calculation shown in the figure the real part 
of the EDAI potential has been used, but a calculation 
with EDADl would give in practice the same result. 

The results displayed in Fig. [3] emphasize the impor- 
tant role of FSI and in particular of the imaginary part 
of the relativistic optical potential, that plays a different 
role in the different approaches. In the rROP the imag- 
inary part is neglected. In the RDWIA it gives an ab- 
sorption that accounts for the flux lost in each channel to- 
wards other channels that are not included in the model. 
In the RGF the imaginary part redistributes the flux in 
all the final-state channels: in each channel it accounts 
for the flux lost towards other inelastic channels and re- 
covered for the inclusive scattering making use of the 
dispersion relations. The results obtained in the differ- 
ent models give an idea of the relevance of these contribu- 
tions. The larger cross sections in the RGF arise from the 
translation to the inclusive strength of the overall effect 
of inelastic channels. We have already noticed, however, 
that while the RDWIA neglects the contributions of some 
channels that can be included in the semi-inclusive reac- 
tion where only the emitted nucleon is detected, the RGF 
is appropriate for the inclusive scattering where only the 
final lepton is detected, and can take into account also 
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FIG. 4. The ratio of NCE/CCQE cross sections as a function 
of calculated in the RDWIA with Ma = 1.03 (solid hne) 
and 1.39 GeV/c^ (dashed hne), and in the RGF with 1.03 
GeV/c^ (dot-dashed-line). The data are from MiniBooNE Q]. 



FIG. 5. The ratio of p/n cross sections as a function of 
calculated in the RDWIA with three different values of As 
and Ma = 1.39 GeV/c^ Results with As = and Ma = 1.03 
GeV/c^ are also shown. 



contributions that are not included in the semi-inclusive 
process. From this point of view the RDWIA can rep- 
resent a lower limit and the RGF an upper limit to the 
calculated NCE cross sections. 

The MiniBooNE Collaboration reported results for the 
flux-averaged differential cross section both in the CCQE 
and in the NCE scattering. In Ref. [l| these results are 
used to extract the NCE/CCQE ratio as a function of 
Q^. This ratio can be useful to compare the results of 
the two measurements. In Fig. 5] we show our results 
for the NCE/CCQE ratio. We note that both NCE and 
CCQE cross sections are per target nucleon, thus there 
are 14/6 times more target nucleons in the numerator 
than in the denominator ^ . The results of the RGF and 
RDWIA with the EDAI potential are similar and within 
the experimental errors. This is in accordance with the 
observation that ratios of cross sections do not depend 
on FSI effects. The RDWIA model, which gives much 
lower predictions for the cross sections than the RGF, can 
produce results for the ratio close to the RGF ones and 
in overall agreement with the data. A more significant 
effect is given by a larger Ma- When Ma = 1.39 GeV/c^ 
the NCE/CCQE ratio is enhanced up to « 30% . This 
means that the axial mass has a different role in the NC 
and in the CC scattering. 

In order to measure the strangeness contribution to 
the axial form factor the MiniBooNE Collaboration used 



the ratio of p/n cross sections for protons above the 
Cherenkov threshold as a function of the reconstructed 
nucleon kinetic energy [l|. In Fig.[5]we show our RDWIA 
results for the p/n ratio. This ratio was proposed as an 
efficient way to measure As [H, [l3l as the distortion ef- 
fects should be largely reduced, but was later given up 
due to the difficulties associated with neutron detection. 
The p/n ratio is very sensitive to the strange-quark con- 
tribution, as the axial-vector strangeness As interferes 
with the isovector contribution to the axial form factor 
gA = 1.26 with one sign in the numerator and with the 
opposite sign in the denominator. We already investi- 
gated in Ref. [l^ the sensitivity of the p/n ratio to As 
as well as to the strange-quark contribution to the vector 
form factors. A large dependence on As was obtained, 
but also the effect of the magnetic strangeness was sig- 
nificant. However, we note that the p/n ratio of Fig. [5] 
is obtained by dividing the flux-integrated cross sections 
with one proton or one nucleon in the final state and, 
moreover, that in the CII2 target there are 8 protons and 
6 neutrons. Thus, it is not straightforward to compare 
the results in Fig. [5] with those of Ref. [l^ . Because of 
the independence of the p/n ratio on FSI, the results with 
the RGF and rROP models are similar to the RDWIA 
ones and are not shown in the figure. The ratio is largely 
enhanced when a negative As is included and suppressed 
when a positive As is considered. Varying the axial mass 
modifies the ratio up to « 10%. 
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III. SUMMARY AND CONCLUSIONS 

In this paper we have compared the predictions of dif- 
ferent relativistic descriptions of FSI for quasielastic NC 
neutrino-nucleus scattering with the MiniBooNE NCE 
data. In the RPWIA FSI arc simply neglected, in the 
rROP they are described retaining only the real part of 
the relativistic energy-dependent optical potential, while 
in the RDWIA and in the RGF the full complex optical 
potential, with its real and imaginary parts, is used to 
account for FSI. 

The RDWIA is based on the same model that was 
widely and successfully applied to the analysis of the ex- 
clusive (e, e'p) knockout reaction, where the final state 
is completely determined. In the exclusive reaction the 
absorptive imaginary part of the optical potential, which 
accounts for the flux lost in the considered elastic channel 
to all inelastic final-state channels, gives a reduction of 
the calculated cross section that is required for a proper 
description of the experimental data. In the RDWIA the 
cross section for the semi-inclusive reaction where only 
the emitted nucleon is detected is obtained from the sum 
of all the integrated exclusive one-nucleon knockout chan- 
nels. In this case, however, the reduction produced in 
each channel by the imaginary part of the optical poten- 
tial, that can be appropriate for the exclusive reaction, 
neglects some final-state channels that can contribute to 
the semi-inclusive reaction. All final-state channels arc 
included in the RGF, where the flux lost in each channel 
is recovered in the other channels just by the imaginary 
part of the optical potential and the total flux is con- 
served. The RGF model is appropriate for the inclusive 
scattering, where only the outgoing lepton is detected, 
and with the use of the same complex optical potential 
it provides a consistent treatment of FSI in the exclusive 
and in the inclusive process. In comparison with data, 
the RGF is able to give a good description of the (e, e') 
experimental cross sections in the QE region and also of 
the recent CCQE MiniBooNE data without the need to 
increase the standard value of the axial mass. The ap- 
plication of the RGF to the semi-inclusive NCE scatter- 
ing can recover important contributions that are omitted 
in the RDWIA, and can give, from the comparison with 
the DWIA results, an indication of their relevance. It can 
also include, however, contributions of channels which are 
present only in the inclusive but not in the semi-inclusive 
reaction. 

The RPWIA, rROP, and RDWIA results generally un- 
derpredict the MiniBooNE NCE data, the RDWIA giv- 
ing the lowest cross section, unless the standard value 
of Ma is significantly enlarged. In contrast, the RGF 
results are in reasonable agreement with the experimen- 
tal NCE cross section without the need to increase the 
standard value of Ma- 

The enhancement of the RGF cross section can be as- 
cribed to the contribution of reaction channels that are 
not included in the other models. It can be due, for in- 
stance, to re-scattering processes of the nucleon in its way 



out of the nucleus, to non-nucleonic A excitations, which 
may arise during nucleon propagation, with or without 
real pion production, as well as to multinucleon pro- 
cesses. These contributions are not included explicitly 
in the model with a microscopic approach, but can be 
recovered, to some extent, in the RGF by the imaginary 
part of the phcnomcnological optical potential. With the 
use of such a phcnomcnological ingredient, however, we 
cannot disentangle the role of different reaction processes 
and explain in detail the origin of the recovered strength. 

If all these contributions can be present in the inclusive 
scattering, the role of multinucleon processes in the NCE 
experimental data is not clear. It is a fact, however, that 
the theoretical analysis of MiniBooNE CCQE and NCE 
data presents a common problem: not only the simple 
RFC, but also other models, based on the lA and includ- 
ing only one-nucleon knockout contributions, require a 
larger value of Ma to reproduce the magnitude of the 
experimental cross sections. The calculations required 
for the theoretical analysis must consider the entire kine- 
matical range of the relevant MiniBooNE neutrino ener- 
gies. Additional complications may arise from the fiux- 
average procedure to evaluate the CCQE and NCE cross 
sections, which implies a convolution of the double dif- 
ferential cross section over the neutrino spectrum. It has 
been argued 0, [2^ that, due to the uncertainties asso- 
ciated with the flux-average procedure, the MiniBooNE 
cross sections can include contributions from different 
kinematic regions, where other reaction mechanisms than 
one-nucleon knockout are known to be dominant. More- 
over, further ambiguities arise for the MiniBooNE NCE 
cross section, which is given in bins where is recon- 
structed from the kinetic energies of the ejected nuclcons. 

Models including other contributions that one-nucleon 
knockout, like our RGF, but also the model of Refs. p?! - 
[29j . where multinucleon components arc explicitly in- 
cluded, are able to describe both the MiniBooNE CCQE 
and NCE data without the need to change the value of 
the axial mass. The two models are different, but they 
seem to go in the same direction. In the RGF, how- 
ever, the enhancement of the cross section cannot be at- 
tributed only to multinucleon processes, since we cannot 
disentangle the role of the various contributions included 
in the phcnomcnological optical potential. 

In order to clarify the content of the enhancement of 
the RGF cross sections compared to those of the lA 
models, a careful evaluation of all nuclear effects and 
of the relevance of multinucleon emission and of some 
non-nucleonic contributions [6l| is required. The com- 
parison with the results of the RMF model, where only 
the purely nucleonic contribution is included, would be 
interesting for a deeper understanding. Processes involv- 
ing two-body currents, whose importance has been dis- 
cussed in Refs. [2^, IH, 113, [l^l , should also be taken into 
account explicitly and consistently in a model to clarify 
the role of multinucleon emission. 

The RGF predictions are also affected by uncertainties 
in the determination of the phcnomcnological optical po- 
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tential. At present, lacking a phenomenological optical 
potential which exactly fullfills the dispersion relations 
in the whole energy region of interest, the RGF predic- 
tion is not univocally determined from the elastic phe- 
nomenology. The differences between the RGF results 
obtained in the present investigation with the EDAI and 
EDAl potentials are visible, alth oug h smaller than for 
the CCQE cross sections in Ref. |52|- These differences 
are produced by the different imaginary part, that is the 
crucial ingredient in both RDWIA and RGF calculations, 
the real part being very similar for different parametriza- 
tions of the optical potential. It is interesting to notice 
that the best predictions in comparison with both CCQE 
and NCE data are given by the EDAI potential, that 
is also able to give the better description of the elastic 
proton-^^C phenomenology. A better determination of 
a phenomenological relativistic optical potential, which 
closely fullfills the dispersion relations, would be anyhow 
desirable and deserves further investigation. 

In this work we have investigated also the role of a pos- 



sible strange-quark contribution As to the axial nucleon 
form factor. The calculated cross sections are almost un- 
affected by As, as the combined effects of As on proton 
and neutron events almost cancel. As a consequence, 
also the ratio of NCE/CCQE cross sections is unaffected 
by As if both proton and neutron events arc included 
in the NCE cross section. Experimental information on 
the strange-quark contribution to the NCE/CCQE ratio 
can be obtained if only proton (or neutron) emission is 
considered. The p/n ratio is very sensitive to the strange- 
quark contribution, but requires the explicit separation 
of the proton and neutron cross sections. 
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